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1 Introduction 



In this paper we generalize the inequahty proved originaUy by Tomboulis and YafFe in SU(2) gauge 
theories [T] to SU (TV) gauge theories with general N and also to a wide range of classical spin systems. 
To make this paper readable for those working on spin systems, we give in this section an elementary 
introduction to studies of quark confinement in lattice gauge theories (LGT) with an emphasis on 
interrelations between concepts in spin systems and those in gauge theories. 

That a specific kind of defect could be responsible for determining a phase structure of a statistical 
system is appreciated as a quite useful idea in wide areas of modern physics. It has a long history, 
possibly dating back to R. Peierls' argument on the Ising model [2]. In as early as 1944 L. Onsager, in 
his famous paper on the exact solution of the two-dimensional square Ising model, calculated what he 
called a 'boundary tension' (the free energy per unit length of a domain wall separating two regions of 
opposite magnetic order) and found that it is zero above and nonzero below the critical temperature 
[3]; hence magnetization is not the only quantity that can characterize the phase structure of the 
system. The history after Onsager clearly tells us the significance of understanding how defects, or 
dislocations, induce the volatility of the order parameter: for instance the seminal work of Kosterlitz 
and Thouless 0] made understanding the infinitely smooth phase transition of the XF-model possible 
by adopting the chemical potential of vortices as an order parameter. 

Such an idea was imported into the studies of SU (N) gauge theories ingeniously by 't Hooft [S] , 
Mack and Petkova [5] and several others [7]. Remember that in spin models, the system is said to be 
in a disordered phase if the two-point correlation function decays exponentially making the correlation 
length finite; otherwise the system is said to be either in an ordered phase or in a Kosterlitz-Thouless- 
type phase. In parallel, a non-Abelian gauge theory is said to be in a confining phase if the expectation 
value of a Wilson loop decays exponentially with the area it spans; otherwise the system is said to be 
either in the Higgs phase or in the Coulomb (or massless) phase. The SU{N) gauge theory with no 
matter field has been believed to be in a confining phase for entire values of coupling constant {SU{3) 
is of special importance as it is supposed to be the true theory of strong interactions in nature, where 
quarks have never been directly observed in experiments). What is mysterious is that a rigorous 
proof of confinement is still missing in spite of a tremendous amount of work dedicated to this issue 
so far. However, according to the scenario(s) pioneered by 't Hooft, Mack, Petkova and others [7], 
the rapid decay of a Wilson loop expectation value might be attributable to a percolation of center 
vortices. It is an object of co-dimension 2 (thus it is a loop in 2+1 spacetime and a closed surface 
in 3+1 spacetime). A rough explanation of their appearance is as follows: in pure SU{N) gauge 
theory, all the fields belong to the adjoint representation of SU{N), so that the actual gauge group 
is SU{N)/Zm rather than SU{N). Iii{SU{N)/ Zm) = Zm means that the systems has a line defect 
associated to each element of Z^^. which is denominated as a 'center vortex', or a "t Hooft loop'. If 
a vortex associated to z G Zn wraps around the Wilson loop, the latter is multiplied by a factor z. 
Then the appearance of infinitely many center vortices piercing the Wilson loop randomly with no 
mutual correlation can efficiently disorder the value of the Wilson loop, resulting in an exponential 
suppression of the expectation value for larger loops followsQ 

Roughly speaking, the formulation of 't Hooft concerns a macroscopically large center vortex 
wrapping around the periodic lattice, ensuring its presence by imposing a twisted boundary condition 
on the lattice. This procedure is essentially tantamount to imposing an anti-periodic condition to 
produce a domain wall in the Ising model. He presented a convincing argument that the behavior 

^It is interesting to note that a similarity of the 't Hooft loop in LGT to the domain wall in spin models gets even 
clearer in the deconfined phase at high temperature. The action of SU (N) LGT possesses a global Zff symmetry, and 
the confinement-deconfinement transition is conventionally interpreted as its spontaneous breaking [8]. The tension 
of an interface separating different Zjy deconfined vacua is calculated perturbatively and numerically from the (dual) 
string tension of the spatial 't Hooft loop [9]. 
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of the free energy of a large vortex in approaching the thermodynamic limit characterizes in which 
phase the system is in; if it vanishes exponentially, then the vortices percolates and the system is in 
the confining phase. On the other hand Mack and Petkova formulated a center vortex contained in a 
torus of finite diameter with a fixed boundary condition on the surface, and the presence of a center 
vortex was ensured by a singular gauge transformation operated on the surface. In order to elucidate 
its intuitive meaning to spin theorists, we would like comment on the concept of 'thickness' of the 
vortex. It is well known that, in the continuum, an infinitely thin center vortex is unphysical in the 
sense that it is associated with an infinite action. For illustration let us consider the XY model on 
a one-dimensional chain of length L. Suppose we fix the angle at one end of the chain to ip and the 
angle at the other end to ip + 9. If the angles of spins change smoothly as much as possible from 
one end toward the other, the energy cost is easily estimated to be {0/L)'^ x L ^ 0(1/ L) for L ^ 1. 
This is in sharp contrast to the situation in the Ising model on the same chain, in which a smooth 
change is impossible, thus leading to the energy cost of 0(1) and making a spontaneous symmetry 
breaking easier to happen. The lesson we learn in this example is that it is generally possible to 
reduce an energy cost associated with a defect by smoothly changing the variables around it; the 
energy cost associated with a domain wall can be reduced if we give it a finite thickness^ What Mack 
and Petkova achieved is to prove an inequality rigorously, whose intuitive interpretation being that 
the area-law decay of the Wilson loop expectation value would follow if the free energy of such a 'thick' 
vortex decreases sufficiently rapidly when its diameter is increased. A lucid exposition of dynamics of 
thick vortices in SU{N) lattice gauge theories (LGT) can be found in ref. [llj. 

As is well known, a fundamental difference between spin systems and gauge theories is that the 
latter has no known local order parameter (such as magnetization in the former) that can character- 
ize the phases of gauge theories. That is why non-local quantities such as Wilson or 't Hooft loops 
have been given a special weight in studies of strong-coupling phenomena such as confinement. As a 
classical reference, we would like to mention ref. 1 12] in which physical relevance of defects generated 
by the twisting procedure, including both center vortices in LGT and domain walls in spin systems, 
and usefulness of using them as a probe for the phase structure in computer simulations, are reviewed 
and discussed from a unified point of view. 

Tomboulis and Yaffc thoroughly investigated SU{2) LGT at finite temperature and rigorously 
proved the absence of confinement at sufficiently high temperature [T]. In their study they derived a 
number of inequalities between observables such as the Wilson loop expectation value, the 't Hooft 
loop expectation value, the electric flux expectation value and the Polyakov loop correlator. Among 
others they gave an upper bound of the Wilson loop expectation value by a specific function of the 
center vortex free energy (as presented in the next section as theorem [T]). It gave a firm foundation to 
't Hooft 's original argument in the continuum [5], that if in approaching the thermodynamic limit the 
free energy of a center vortex that encircles two of the four periodic directions of the lattice vanishes 
exponentially w.r.t. the cross section of the lattice perpendicular to the vortex, then the area law 
behavior of the Wilson loop expectation value would follow. Thus it sheds light on dynamics of the 
center vortices in a somewhat different manner from the Mack-Petkova inequality. In this paper we 
call it the Tomboulis- Yaffe (TY) inequality throughout this paper. 

The purpose of this paper is to present a generalization of the TY inequality to SU{N) LGT for 
arbitrary N and to general classical spin models. Our result gives a rigorous relation between the 
effect of twisted boundary conditions and the correlation function (Wilson loop) in spin models (in 
LGT) respectively!! Among spin models, the SU{N) x SU{N) principal chiral model (PCM) is of 

^ Dobrushin and Shlosman elevated this idea to a rigorous proof of the absence of magnetic order in two-dimensional 
ferromagnets with a continuous symmetry 1101 . 

•^Historically, changing of boundary conditions has been utilized in studies of Anderson localization as a method for 
estimating the broadening of the wave function 1131 . More recently it was utilized in the lattice QCD calculation I14| 
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particular interest for researches of gauge theory, since it bears a number of similarities to SU{N) 
gauge theories and serves as a good testing ground for techniques in gauge theories [13 [HI [T71 [TH] . 
The action of SU{N) PCM is given by 

d 

S = (3^^ReTT{U{x)U''{x + fi)}, U € SU{N), x € Z'^, (1) 

X /J — 1 

where fl denotes a unit vector in x^-dircction. It is quite straightforward to extend the original TY 
inequality for SU{N) LGT to SU{N) PCM, using a natural correspondence (site <-> link, link ^ 
plaquette, . . . ) and indeed the TY inequality for SU (2) PCM has already appeared in the literature 
[l9l [20] , On the other hand, however, it is technically nontrivial how to extend it to other more 
general spin models. Let us take G2 PCM as an example. Since G2 is an exceptional group with 
trivial center, we can no longer use a twist by a center of the gauge group, which gives rise to a 
technical difficulty. Furthermore the use of center twist for PCM is not physically motivated; in the 
case of SU{N) gauge theory, the use of center element is mandatory, but in PCM we can use any 
other element of the symmetry group for twist. Thus the generality of our formulation, that does not 
rely on the center of the symmetry group at all, seems to be a fundamental progress|j 

This paper is organized as follows. In section 2 we will recapitulate the TY inequality for SU{2) 
and then prove its generalization to SU{N). We will use the two-dimensional SU (iV) LGT to illustrate 
our result. In section 3 we will prove a generalization of the inequality to general classical spin systems. 
We will use the one-dimensional PCM and the two-dimensional Ising models on square and triangular 
lattices to illustrate the proved inequality. Especially, in section lXSl we derive a rigorous upper bound 
of the off-axis correlation length in the triangular Ising model, whose exact expression is still unknown, 
and conjecture that it is indeed the exact one. In section [XHl the strong coupling expansion technique 
is employed to shed light on the implication of our theorem, as well as to test the conjecture. Section 
4 is devoted to the conclusion. 



2 TY inequality in LGT 



2.1 N = 2 

Let us recapitulate the TY inequality for SU{2) [T]. A is a c?-dimensional hypercubic lattice of 
length Lf^ = 1, . . . with periodic boundary condition and V, called "vortex", is a stacked set 
of plaquettes winding around the lattice A in 0? — 2 periodic directions!! We assume the directions 
unwrapped by V to be x^^ and x'^ (/i 7^ ly). The ordinary and the "twisted" partition functions are 
given by 



/3 



2 ETrC/p), 

pCA 



UdU 



b exp 



2 



ETr(-^p) 

pCV 



pCA\V 



(2) 



(3) 



where dU is the normalized Haar measure of SU{2) and Up = Ux,^Ux+fj.,iyUl_^_^ ^U^ ^, is a plaquette 
variable. It is important that local redefinition of variables U —U can move the locations of twisted 
plaquettes but cannot remove the twist from A entirely. 



to study charmonium properties in deconfinement pliase. 

^ As an aside we note that the incquahty of Mack and Petkova for SU (N) LGT was generalized to SU (N) PCM by 
Borisenko and Skala |21) . 

forms a closed loop when d = 3 and a closed surface (2-torus) when d = 4, on the dual lattice. See fig^ 
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Next, consider a rectangle C lying in a x^-x'^ plane with Ac the area enclosed by C, and let W{C) 

the Wilson loop in the fundamental representation associated with C, namely W{C) = "^Tr JJ^ Ub- 

bee 

Then the following inequality holds [TJ [52] : 
Theorem 1. 



1 /. Z 



(-) 



WC))<2|-(^l-^j| (4) 

where (...) is the expectation value w.r.t. the measure of Z\. 

The site-reflection positivity of the Wilson action [23] plays an essential role in the proof. (As is well 
known, the Wilson action is among those actions for which the link-reflection positivity is also satisfied 
fM\ but it is not a matter of interest here.) Indeed ([U) can be proved with any one-plaquette action, 
since they are site- reflection positive (although not necessarily link- reflection positive, of course). 

An important implication of ^ is that the area-law decay of {W{C)) would follow — V-^A ^ 
Q-pL^L^ /or some constant p > in the thermodynamic limii^ This is a famous criterion of confine- 
ment originally proposed by 't Hooft |2] and is also numerically supported [55] [IS] [57] . Moreover such 
a behavior of .Z^ V^A has been verified explicitly by Miinster using the convergent strong-coupling 
cluster expansion technique [5Sj. See page [TU] for more discussion on this point. 

Theorem [T] was utilized in a recent attempt at a rigorous proof of confinement [29| with related 
discussions [30l [311. 



2.2 General N 

The authors of ref.[T] state without explicit construction that their result is extendable to any other 
gauge group with nontrivial center. Since the mentioned extension does not seem to be so trivial and, 
to the author's best knowledge an explicit formula for general N is not found in the literature, we 
think it valuable to present the extension of ^ from SU (2) to SU {N) together with its proof. 

Let us give a formulation of vortices in SU{N) LGT and prove their properties before presenting 
TY inequality in SU{N) LGT. The ordinary and the "twisted" partition functions are respectively 
given by 

Za = fl[dUbexpl^Yl ReTrC/p , (5) 

h \ pCA J 



Z^'^JUdUbexJ^ 



^ Re Tr(z'= C/p) + ^ ReTrC/p 

pCV pCA\V 



(6) 



z = exp(— j, A; = l,2,...,iV-l(modiV). (7) 

Hereafter (...) represents the expectation value with the measure ([5]). The vortex creation operators 
{C'[*''l} and the electric fiux creation operators {J^I™!} are defined by 

[k] 

(OW[V])^^, (8) 



1 ^-1 

(^[™l[v])^_^^"'^(OW[V]). (9) 



k=0 



°We neglected the entropy factor for simplicity. 
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N-1 



Thus we have (C'['=1[V]) = ^ [V]). The expUcit form of O is given by 

m=0 

[V] = exp(-^ P° Tr (z'^'C/p) - Re Tr Up]^ . 

pCV 

Lemma 1. //V, V", V", . . . are homologou^, we have 

(C)[fel[V]C)['^'l[V']) = (C'[''+'^'1[V]), 

m=0 

(^W[V]^[™l[V']) = (^['l[V])4fj, 



(10) 

(11) 

(12) 

(13) 
(14) 



where S^^J = 1 if I = m (mod N) and ~ otherwise. 

Proof. (fTT|) can be derived by iterating the redefinition of variables U ^ z'' U to bring V" to V. 
Relations (HHl-dlll) follow from (O and □ 

P2p - (|14p imply that {J-^"^^}m can be seen as projection operators [SI 111], 

Consider a (d — l)-dimensional hyperplane tt defined by x'^ = m with m G Z fixed! Define the 
reflection operator 6 w.r.t. tt by 0[F({[/h})] = F{{Ug[ij]}) where F is an arbitrary observable (that 
is, a map from configurations on A to C). The reflection 9[b] of a link b is also defined by the same 
notation where locations of b and 0[b], are defined to be symmetrical about tt. 



Lemma 2. With V* = 0[V] 



jrM[V]j = J^["1[V^], 

< (0[''1[V]) < 1, 
< (Jc-H[V]) < 1. 



(15) 
(16) 

(17) 
(18) 



Proof. (fTS|) is obvious from the fact that the orientation of plaquettes are reversed by reflection. ([T5|l 
yields 



1 



JV-l 

k=0 
N-1 

fe=0 



which proves (jl6p . Next, using the Schwarz inequahty < (FOPy/^ and (HIl), (US]) we find 



(o['=i[v]) < (oW[v]0b['=i[v] 



,1/2 



(19) 

(20) 
(21) 

(22) 
(23) 



^Plural vortices are called homologous if and only if they wind around the same periodic directions of A. 
* In this paper we never use hyperplanes defined by x'' = m + ^; that is, we never use link-reflections. 
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which proves the second inequahty in pT|) while the first one is trivial. Since V and are homologous 
we can apply (jl3p to obtain 



(^H [V]) = (^H [V]^['"l [V^]) (24) 

= (jC-HiViefjc-I^liV]]) > 0. (25) 

and ini) prove (HH). (These simple proofs of and (HH) seem to be new.) □ 

The vortex free energy i^i'^' and the electric flux free energy F^f are defined by e^^-" = (0['^1[V]) 
and e~"^<i' ' = (JF[™1[V]); respectively. 

Let N{R) e {0, 1, . . . , - 1} denote the A^-ality of an irreducible representation of SU{N) 

whose dimension is dji. Take a rectangle C lying in a x^-x"^ plane and let Ac the area enclosed by 

C. For the normalized Wilson loop in the representation R, VFfl(C) = — — xb.( Y\_ ^'>) ' have 



bee 



Theorem 2 (TY inequality for SU{N) LGT). 

\{Wr{C))\ < (^[JV(i?.)][v])Ac/L,L„ ^ 1^ _ (^™[V])} 

/n addition, if N{R) ^ we have 

|(W^i?(C))|<2{l-(.F["l[V])}^"^'^''''" 
f 1 



(26) 

(27) 
(28) 



k=0 



Proof. Although the argument below parallels that of ref. [22 for SU{2), we describe the proof in full 
detail for readers' convenience. Suppose V, V are stacked set of plaquettes wrapping around d ~ 2 
periodic directions of A and V is linking once with C while V" is not. (See figlDfor a 3-dimcnsional 
illustration of the setting.) 




Figure 1: Locations of V, V" and C. 
Let us rewrite the expectation value of Wji{C) = -i—Xb ( Ub) as follows: 



bee 



{Wr{C)) = ((1 - .Fl"l[V])M/fl(C)) + {T^"^[V]Wh{C)) 



= ((1 - mWniC)) + ^ E MWBiC)) 



(29) 
(30) 



fc=0 



'Af-ality is the number (mod A'^) of boxes in the Young tableau of R. 



The presence of the second term is not desirable from the viewpoint of obtaining a meaningful upper 
bound of {Wr{C)), so let us perform redefinitions of variables U z^U to bring V to V", which 
causes the change 

{0^^^[V]Wr{C)) ^ z^^^^^^{0^%V']Wr{C)). (31) 

This is because in the course of bringing V to V' we must change one of the link variables on C . (The 
sign of exponent depends on the orientation of C.) Thus pop becomes 

{Wr{C)) = ((1 - ^["l[V])Ty,^(C)) + ^ E ^^''^"''\0^''\V]Wr{C)) (32) 

= {{1 - T^"^[V])WRiC)) + (^l±^(^)l[V']VKi^(C)). (33) 
Our next step is expressed in figH] schematically in which a black square represents the oper- 

Figure 2: The Schwarz inequality enables us to double the size of the rectangle. This figure is 
essentially borrowed from ref. |22j . 

ator (1 — jrM[V]). Labeling four vertices as A, B,C, D, setting the hyperplane tt so that it is 
perpendicular to the rectangle and contains the edge CD, and applying the Schwarz inequality 
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\{FdG)\ < {FdPy/'^iGeGy/'^ we obtain 

\ii^-^^"^mWn\DA+AB+BC+CD)\ (34) 

= ^|((1--^'°1[V])X«|^^^^^^^^^^^)| (35) 



a, [3 



^ E ((1 - - ^'°'[V']){x«L^+^^^^c}"/3{x«L^.+^.^«+^.c:}"/3>''' 

/ j \l/2 

^{Xfl|c£,}/3a{Xi?|c£,}/3a/ (37) 

= {{l-^^°HV])Wn\^s^sss+se^e^^e^)'/'. (39) 

In going from ((57)) to ([55]) we used The length of the rectangle doubled. 

Let Lx li, the size of the original rectangle and assume that = 2^1^ and = 2'/^ for some 
p,q ^ INlij. Repeating the operation above for sufficiently many times in both x^- and x"- directions, 
due to the periodic boundary conditions the rectangle finally vanishes away, yielding 



{{1-T^°\V])Wr{C)) <[l- (^[°1[V])} (40) 



^jr[±w(fl)][V']M/fl(C)) < (jc-[±w(«)][v'])Ac/-L^,L. can be shown in a similar way, hence ^ is proved. 
dUl) is a consequence of (HH) and [V]) < 1 - [V]). □ 

The message of (|28p is that the exponential decay of the vortex free energy, i.e. (O^'^'l [V]) = c^^" ' = 
1 — 0(e~''^'^^") for every k, is a sufficient condition for the area law of the Wilson loop to hold. The 
area law docs not hold, or is at least difhcult to prove, if not all of the (C')'s converge to 1. 

Several comments are in order. Firstly, suppose that the action in ([2]) is in the adjoint representa- 
tion. Then (C'[''1[V]) = 1 follows for any k, hence making {Wr{C)) with N{R) ^ vanish identically 
for arbitrary finite volume (sec (pS]) ). This is to be anticipated; since the adjoint Wilson action is 
invariant under the local transformation U — > zU (z G Zn), which cannot break spontaneously ac- 
cording to the Elitzur's theorem, and since Wii{C) with N{R) ^ transforms nontrivially under this 
transformation, its expectation value must vanish. 

Secondly, there are a lot more varieties of inequalities available other than (|26p . Assume A'^ = 6 
for instance. From ^ we have ((JpW[V] + T^^^[V]){T^^^[V] + T^'^^[V'])) = {T^^^[V] + T^^^[V]) and 
( (J^[3] [V] + W [V] + ^[51 [V] ) (^[3] [ v] + jr[4] [V] + ^[5] [ V'] ) ) = (.^[3] [V] + [V] + .FI^I [V] ) , hence by 
modifying the above proof one can straightforwardly show 

However (j42|) and all of its cousins are weaker than ((26|) with N = 6, which can be understood by the 
elementary inequality (^- Xi)" < X]i(^i)" for < a < 1 and < Xi. 



^"This condition was also present in the original TY inequality It is not a severe restriction, however, as long as 
we believe the asymptotic behavior of observables to be independent of the way we take L^, to infinity. 
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Note that one cannot derive the area law from (pS)) when N{R) = 0, as can be seen from 
r.h.s. of ^ 

The above hiiphes that the "ghions" of SU{N) can screen particles of zero iV-ality. 

Thirdly, theorem [5] is correct even after a matter field whose 7V-ality is zero is introduced into 
the theory, since the matter-gauge couphng ^l.^p^Dr[Ux,fj,]^x preserves reflection positivity and is 
insensitive to the change of variables U ^ zU. 

Finally we remark on the utility of strong-coupling cluster expansion techniques. (Similar discus- 
sion will be presented in section [3.61 ) As already mentioned, the exponential suppression of vortex 
free energy -log(e'[V]) w c'P^'^^" has been verified by Miinster [H] for SU{2) LGT and for suf- 
ficiently strong coupling. Especially he showed to all orders of strong-coupling expansion that the 
constant p appearing in the vortex free energy ('t Hooft's string tension) is equal to the conventional 
Wilson's string tension. His proof hinges on the observation that both the calculation of Wilson loop 
expectation value and that of vortex free energy reduce, at sufficiently strong coupling, to the problem 
of fluctuating random surfaces. It is understood without difficulty that the methods he employed can 
be readily used for SU{N) LGT to show - log(0['^l [V]) « e'P^f^^" for every fc 7^ (hence proving the 
area law) . This time p is equal to the fundamental string tension (since the gauge action ([5]) is in the 
fundamental representation). 

* * * * 

Let us then turn to LGT with matter field of non-zero N-ality; the relations pT|) - (fH|) are no 
longer valid. If the matter field has A'^-ality m and the greatest common divisor of N and m is s , the 
subgroup Zs C Zn is a symmetry of the theory. It is thus straightforward to prove the following 

Theorem 3. If N{R) ^ (mod s), we have 

\{Wn{C))\<2{\-{f^^\V])]" '-=2|l--^^^— I . (44) 



>(^[o][v])Ac/L.L„>/ M — ^1 as L^,L,^oo. (43) 



Although the development so far has been for SU{N) LGT, the inequalities evidently apply to 
U{N) LGT since the center of U{N) is U{1) which contains all of Z2, Z3, Z4, Z5, .... Let us focus 
on C/(l) for simplicity and define the twisted partition function as 



Za{0) ^ Jl[dUb cxp /3 



^Rc (e''^ Up) + Rc Up 

pCV pCA\V 



We state below the counterpart of theorem [H The proof is straightforward. 
Theorem 4. For the Wilson loop of U{\)-charge q E \ {0}, we have 



I(w^,(c))|<2{i-(.f[%[v])} 



Ac/L,,L,, 



de 
"27 



ZA{e) 



with 



27r 



de Zjy{e) 

"2^ Za 



(45) 



(46) 



(47) 



Finally we point out that theorem [2] can be proved even if we add a Wilson loop with zero N- 
ality, of size 2 x 2 or 2 x 1 to the action. It is simply because the site-reflection positivity is kept and 
the algebras of twists are still well defined. 
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2.3 Demonstration in 2D SU{N) LGT 

Let us explicitly verify the proved inequality in solvable two-dimensional SU{N) LGT. In two dimen- 
sion, the twist is introduced on just one plaquette. Let the size of the lattice Li x L2 and impose 
periodic boundary conditions in both directions. First we expand the exponentiated one-plaquette 
action, e^'^p, into characters of irreducible unitary representations of SU{N): 

e-S'^iU) =Y,drFrXr{U) (48) 

r 

where denotes the dimension of a representation r and the reality of Sp implies Fr — F- (ovcrline 
represents complex conjugation). Then a straightforward calculation using formulae 

dUXr{VU)Xr'{U^W) = j-5rr'Xr{VW), (49) 

j dUXriVUWU^) = l-Xr{V)Xr{W), (50) 

yieldf0 

Zk = jY{dUb\{{Y. drFr-XriUp) J (51) 

beA P \ r J 

= Y^{Frf^^\ (52) 

r 

On the other hand, introducing a twist = exp(2Trik / N) into arbitrary one plaquette on A gives the 
twisted partition function 

4^=1 =^(F,)^i^-z'^^W. (53) 



Af-l 



Using the identity ^ z^^i'^) = iV(5o_Ar(r) we easily obtain 



fc=0 



N 

k=Q 



(55) 



i + E(^'-)''^^^ 

Fr 

where T implies the trivial representation and we defined = -= — . Note that Cr = Cf. Since 

Fx 

\Fr\ = J dU e-^'^^^^UU)\ < J rfJ/c-^p(^) =Ft, (56) 



^^1(52} differs from that obtained in ref . |32) because tfiey impose free boundary conditions. 
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we have \cr\ < 1, while < Cr can be shown for a wide class of gauge actions including the Wilson 
action. 

From above considerations we obtain 

r N-l ~| ^c/LiL2 

r.h.s. of (EH)] = 2 1 - — 5^ (O W ) (57) 
I fc=o J 

-> 2(cr/)^^ as Lii2 ^ oo. (58) 

Here Cr' is defined as the largest value among {cr \ N{r) =/= 0}. In order to determine r' we need an 
explicit form of the action Sp. 

Let us turn to the l.h.s. of ((27|) . i.e. the Wilson loop expectation value. We borrow the result of 
ref . [32j which in our notation reads 

(W^«(C)) = (c^)^-. (59) 
If N{R) ^ 0, we obviously have < c^', hence the inequality (|27p holds for sure. 

3 Extension of inequalities to spin systems 

Main result of this section is theorem [5] on pagcfTTl which is a generalization of theorem [2] to general 
spin systems. Before that, we need some preliminary analyses. 

3.1 Basic formulation 

Let us formulate a twisting procedure in spin systems obeying ref. |12] . Consider a statistical system 
with nearest-neighbor interactions whose partition function is given by 



j W dcf>^ cxp I ^ ^ A{c^y,q^y+f,) j , (60) 



where A is a d-dimensional hypercubic lattice with periodic boundary conditions, /i is a unit vector 

d ' 

in the /i-direction and is a sum over all links in A. The real-valued symmetric function A{ , ) 

yeAfj.=l 

dictates the interaction between nearest sites (and possibly includes self interactions on each site). It 
can be shown by standard arguments that site-reflection positivity is automatically satisfied for any 
nearest neighbor interaction (sec p. 33 of ref. [23]) while link-reflection, not needed in the following, is 
often violated. Hereafter (...) represents the expectation value with the measure ([60|) . Let us assume 
that the system is invariant under a global transformation (j) gcp for any element g of a global 
symmetry group 

A(<^,0')-^(50,.90'), 9^G. (61) 

We assume that G is compact. 

A twist for a link is defined as the change of interaction from A{(j), (j)') to A{(j), gcj)'). An important 
difference from the twist in LGT is that g need not belong to the center of G. G may or may not 
have a nontrivial center and that is not important for us. 

Next, let us take a stacked set of links, V, which winds around A in d — 1 periodic directions (V is 
a closed loop when d = 2 and a closed surface when d = 3 on the dual lattice, see fig. E]). Hereafter 
such V is called a wall in distinction from a (center) vortex. 



^■^Note that G need not be the maximal symmetry group of the system. The development in this section still holds 
if we take as G an arbitrary subgroup of the maximal symmetry group. 
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Figure 3: V in two- and three- dimensions. 



The twisted partition function associated to g G G is given by 

d 



■(62) 

The symbol (5[V] y.^i is defined to be = 1 if the hnk {y, y+fi) is contained in V and = otherwise. It is not 
difficult to see that V cannot be removed from A by local redefinition of variables g4>. Let G denote 
the set of irreducible unitary representations of G. The wall creation operators {0(5) [V]| g S G} and 
their duals {.?^ii:((7)[V]| g E G, R E G} are defined as follows: 

(0(5)[V])^^^|M, (63) 

d 

^ O(.g)[V]=exp(^^<5[V],,^[-^((/.,,0,+^)+A(0„50,+^)]), (64) 
J-fl(g)[V] = (dimi?) / dxO{gxmxR{x). (65) 

JG 

Lemma 3. If the walls V, V" are homologous^^, we have 

{Oig)[V] ■ Oig')[V']) = {0{gg'm) = {Oig'g)[V]), (66) 
{Tn{g)[V] ■ J^R'ig'm) = SRn'{:Fn{gg'm) ^ SnR'iJ^Rig'gM), (67) 

{Oi9m)^J2(^MV]), (68) 
Red 

l=Y.{Tnil)[V]). (69) 

Red 

i66]) . ( [6'7| j imply that {0{g)[V]) and {J-ii{g)[V]) are class functions on G. 

Proof. (j66p is trivial, since the relative position of walls can be reversed owing to the periodic boundary 



^^Plural walls are called homologous if and only if they wind around the same periodic directions of A. 
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condition (see fig. [3]). (p7|) can be shown by exploiting the invariance of the Haar measure: 
(^i?(5)^i?'(ff')) = (dim i?) (dim i?') ( f dx 0{gx)xR{x) [ dy 0{g'y)xR'{y) 



\JG 

= (dimi?)(dimi?') [ dx ( dy {0{gxg'y))xR{x)xR'{y) 

JG JG 

= (dimi?)(dimi?') / dx{0{x)) f dyxR{i99T'xy-')xR'(y) 

JG JG 



SRR,{dmiR) / dx {0{x))xR{{gg')-'x) 

JG 

Srr' {J^Rigg'))- 



(70) 

(71) 

(72) 

(73) 
(74) 



FinaUy, ((55)) is a consequence of the Peter- Weyl theorem according to which any / e L?{G) can 
be represented as 

f{g) = (dim^) / f(^9)XR{x). (75) 



ReG 



triviaUy follows from 



□ 




9 g 



Figure 4: Relative position of V and V' can be reversed by using the periodicity of A. 

Note that (|66p cannot be proved in general when other operators are inserted, that is, 

(0(g) [V] • 0{g')[V] + {0{gg'W\ . . .) + (0(g'.g)[V] . . . ), (76) 

in general, where . . . denote additional insertions. It is because the proof of (j66|) involves a sequence of 
changes of variables. Thus (j66p could be proved if inserted operators are invariant under the changes 
of variables. 

Take an arbitrary (d — l)-dimensional hyperplane tt defined by x'' = m, m ^ X with m and /x 
fixed. Denote by 9 the reflection about tt. 

Lemma 4. // the wall V seen on the dual lattice is also perpendicular to the x^-axis (see the left of 
fig. we have 



0[Oig)[V]\ ^ 0{g-')[V% 
{Oig)[V]) = {Oig-')[V]), 



{TR{g)[V]) = {Tnig-^m) and (^fi(l)[V]) £ R, 

ol^Rigm] ^:Fn{g-')[V% 



(77) 
(78) 
(79) 
(80) 



where = 0[V] and 1 denotes the unit element of G. 
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TT 



y 


z 


• — • 

y + A 


• — • 

z + p. 




> 



Figure 5: Insertion of a reflection plane tt. 

Proof. From the right of fig. [HI we observe 

e A{(i)y,gcj)y+f,) = A{(l)^+(,,gcj)^) 

The twist g changed to g~^ , hence ([77| is proved. (j78p immediately follows from ([73 
Using ([75]) we can show 

dx {0{gx)[V])xB{x) ^ I dx {0{{gxr^m)xB{x) 

G JG 

dx {0{x-^g-^)[V])xR{x) 
dx {0{xg-^m)xR{x-^) 



G 



G 



dx {0{g-'x)[V])xR{x), 



(81) 
(82) 



(83) 
(84) 
(85) 
(86) 



therefore ((7^ is proved. In the last step we used the fact that i? is a unitary representation. (|80p is 

obvious from ([7^ . □ 

Though it seems hard to find more properties on general grounds, further nontrivial result can be 
obtained if we exploit the site-reflection positivity of the measure of (pUj) . 

Lemma 5. 

(0(ff)[V]), {Tn{l)[V]) e[0,l], (87) 

|(^fl(5)[V])| < (^fl(l)[V]). (88) 

Proof. Let = 0[V]. With the aid oi^,^ and the Schwarz inequality < {F9Fy/'^ we get 

{0{g)[V]) < {Oigme[o{g)[V]\y/' (89) 

= {0{gmO{g-')[V'])'/' = 1. (90) 

Next, using and ([80|) yields 

{:FRii)[v])^{Wm-W)[v']) (91) 

= (J-fl(l)[V]-0fj-fl(l)[V]l)>O. (92) 
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(PI) combined with ^ yields dST]). 

Finally, to show ^ we use ^ and the Schwarz inequality \{FeG)\ < {FeF)^/'^ {060)^/'^ as 
follows: letting V' denote a wall homologous to V, we have 



I(-^«(5)[V])| = |(^h(,9)[V]^«(1)[V'])| (93) 

< (^«(5)[V]^K(ff-^)[V^]>^/'(^fl(l)[V']^H(l)[V"^])'/' (94) 

= (^«(l)[V])i/2(^^Xl)[V'])^/' (95) 

= (■^i?(l)[V]). (96) 

□ 



3.2 TY inequality in spin systems 

A natural counterpart in spin systems of Wilson loops in LGT is a two-point correlation function 
(py) as explained in the Introduction. V will decay exponentially with a mass gap (in symmetric 
phase) while decay algebraically without a mass gap (in a spontaneous symmetry breaking phase or 
Kosterlitz-Thouless-type phase). If {0{g)) defined above converges to 1 in the thermodynamic limit, 
it follows that arbitrarily huge domain walls grow with little cost and eventually drive the system to 
the disordered phase with a mass gap (see fig. [5]). {0{g)) 1 (|A| oo) can also be regarded as 
a sign of insensitivity of the system to boimdary conditions, which indicates the absence of massless 
particles. 

If (as in Ising-like models) we assume that the intersection of a correlation line with a wall changes 
the sign of F, a small closed wall gives no contribution ((—1)^ ~ 1) while a huge wall can give (— 1). 
If each link on the correlation line of total length L is assumed to intersect with a wall independently 
with probability p, we obtain 




Figure 6: Disorder being caused by walls. 

Our goal in this section is to elevate the above heuristic relation of disorder and a mass gap to a 
mathematically rigorous inequality. 

Suppose that the explicit form of the correlation function F is given by 

F^(a;; n) = ^ fa{<Px)fa{(l>x+nii), n e N, (98) 

a 

where / is an arbitrary function from the order-parameter space to C endowed with generic indices 
{a}. Here fa is meant to specify, for example, the a-th component of a vector spin in 0(A^)-like 
models, or the a-th matrix element of a matrix spin in PCM-like models. Note that (Fp(a;;n)) is 
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independent of x due to the translational invariance of the system. An important requirement on 
is its invariance under G : 



^fa{(f'x)fa{(l)x+nfL) = ^ f aigcj^x) f a{g(l>x+np.) , 9 
a a. 

Another requirement, which is truly indispensable for all the following development, is 

/ dgU9(t>)=0 for 

JG 

If ()100p were not satisfied, we should replace fa{4>) by /4((/>) = fa{4') ^ / dg fa{g4>)- 

JG 



(99) 



(100) 



Theorem 5. Assume 
of A in the -directioni^ Then we have 



\), ilOO\) and the existence o/ fc G N such that = 2^n, with the extent 



< 2 



(r,(x;0))2 



dgiOigM) 



(101) 



where V is a wall perpendicular to the x^ -direction 

Proof. We first decompose the correlation function into two parts: 

{T^{x;n)) ^ {T^{x;n){l - Triim)) + {T^{x;n)TTiim), (102) 

where T denotes the trivial representation, i.e. J^t(1)[V] = / dgO{g)[V]. Our basic idea here is to 

JG 

apply the Schwarz inequality 

\{FeG)\ < {FeFY/^{GeGY/^ (los) 

to each term of ()102p . The procedure afterward is represented graphically in fig. [7]where (1— JFr(l)[V]) 
is indicated by red segments and a blue line is drawn to guide the eye. Let y = x + nfl and z = x+2np,. 
Consider a reflection 9 about n (a hyperplane perpendicular to fi and lying a.t y ). Using ([^ . ([501) , 
([M]) and pll^l) we get 



\{r,{x;n){l-TTil)[V]))\ 



(rM(2^;0)) 



< 



1 



(rM(2^;0)) 



\J2ifMfMii-TTiim)) 



1/2 



U{(t>v)0 fa{4>v] 



< 



■ E - ^t(1)[V])/.(0.)(1 - .fT(l)[V'^]))'^"(/„(0,)/.(<^,)) 

a. 

E )x:(^(i ~ .ft(i)[v]))'^V E ^(0;)/„(0,)) 



(rM(a^;0)) 
1 

(rM(a:;0)) 



(r^(x;2n)(l-.Fr(l)[V])) 
(r;.(a:;0)) 



1/2 



(104) 
(105) 

1/2 

(106) 
(107) 

(108) 
(109) 



It does not seem to be very restrictive; in general wc believe in the existence of the limit L^^ — + oo for physical 
observables independent of the way we take —> oo. 

^^0 < (Tfi(n)) can be proved by the site-reflection positivity if n is even. 
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Figure 7: Schwarz inequality enables us to double the distance of points in the correlation function. 



Iterating this procedure for (fc — 1) times yields (note ~ 2^n) 



|(r^(n)(l-.FT(l)[V]))| 



< 



(r^(£^/2)(l-.FT(l)[V])) 1 

(r,(o)) / 



(110) 



Finally, let us define the reflection 9 w.r.t. the hyperplane tt which runs through x and x+{L^/2)fj, 
(thus 9[x] = X and 9[x + (Lp/2)/i] = x + (i^/2)/i, see figlHl)- Then we obtain 



< 



< 



{r^{x;L^/2){l~TT{im)) 

(r^(x;L^/2)0[r^(x;L^/2)]}V^((l-.FT(l)[V])(l-.FT(l)[V^]))^/^ 

(r^(a^;0)) 



(rp(x;0)) 



(rp(^;0))^ 



{1-(J-t(1)[V])} 



1/2 



thus 



< 



{T,{x;n)il-TT{im)} ^ f (r^(x;0)2) ^-Z^- 



(r;.(x;0)) 



{i-(^T(i)[v])y 



(111) 

(112) 
(113) 

(114) 
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V 



e 



Figure 8: is a reflection about tt that bisects A. Here A is shown as a one-dimensional chain for 
simpHcity. 

In deriving pl^ we used (r^(.T; L^/2)6'[r^(.T; i:^/2)]) < (r^(x;0)2). This can be shown as follows: 



{T^{x-Lj2)e[V^{x;L,j2)]) 

X! Uo,{<Px)fa{(px+(Lj2)i,) ■ .fp{4'a;)f0i(l}x+{L^/2)il) 



X] {fai(l^x)fa{<Px+{L^/2)t,) //3 (0x )/,3 ((/'a:+(L^/2) A ) ) 



(115) 

(116) 
(117) 



1/2 



{<Px+{L^/2)p.).fl3{4'x+{L^/2)fL) ■ 0\fa{<Px+{L^/2)fi)fl3{4'x+{L^ 



< (Ei/"('^-)i'i//'('^-)i')'^'(Ei/"('^-+(i./2)M)ri//3('/'.+(L./2)A)i'''^' 

a, (3 a, 13 

= {T,{x; OfY'^{T,{x + (L^/2); = (r,(x; 0)^). 

(ru(a;;n)J^T(l)[V]) 
Next we have to estimate '"^ ' ' ^ j/ 



(118) 
(119) 

(120) 
(121) 



. The outline is similar to the previous case, but 
additional intricacies occur. 

0< (r^(x;n)^T(l)[V]) (122) 
= E / d9{Uct^.,)U^)0{g)[V]) (123) 

Let us consider another wall V" that is parallel to V and bookends x with V (see figE])- Moving V to 
V passing over (ji^ and using (|100|) . we get 

(124) 
(125) 



^ / dg{Ugcl^,)Ucj,y)0{g)[V]) 
E / dg(^Ugct^,)JM(0{g)[V]~TT{l)[V] 
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This step (from (|124p to (|125p ) is the most nontrivial operation in this proof. Using (|103p w.r.t. the 




1/2 



Figure 9: By redefinition of variables. V is moved to V', passing over the site x. 
hyperplane which runs through y, with z = 0[x\, yields 

(126) 
(127) 

= E/ dg{^fM.)lM7)[l-0{g)[V]TT[l)[V']-0[g-')[V'^^^ 

< (r^(2/;0))i/2 / dg(r^(x;2n)(l-0(g)[V']^T(l)[V'V0(5"')[V'VT(l)[V']+^T(l)[V'] 



1/2 



(128) 

1/2 



< (r^(x;0))i/2j / dg^r^(^;2n)(l-0(5)[V']^T(l)[V'V0(.9"')[V'VT(l)[V']+.FT(l)[V'])) 



(129) 

1/2 



(r^(x;0))i/2r(r^(^;2n)) - {r^{x;2n)TTil)[V']TT{l)[V'']) - (r^(x; 2n).FT(l)[V' Vt(1)[V']) 



(r^(x;2n).FT(l)[V']) 



1/2 



(r^(x; 0))i/2(r^(x; 2n)(l - ^T(l)[V']))l/^ 



(rM(-^;0)) 

Combining this result with (|114p yields 



(r^(x;n).FT(l)[V]) ^ / jry.T; 2n)(l - .Ft(1)[V'])) 



(rp(a;;0)) 



1/2 



(r;.(x;0)) 
(dm) and (fT34l) lead to pm|) . 
Here are a few comments: 



(r^(x;n)^T(i)[v]) < r (r^(.T;0)^) ^ "/"^^ 



{l-{TTil)[V'])} 



/111 n/Lj, 



(130) 

) 

(131) 
(132) 

(133) 

(134) 
□ 
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• {0{g)) is a class function on G. It implies that when G is a compact connected Lie group we 
can rewrite the integral in the r.h.s. of (|10ip as an integral over the maximal torus of G with a 
proper Jacobian, using the Weyl's integration formula 

• According to pOip . the exponential decay of the wall free energy, i.e. 

(0{g)[V]) = 1-0(0- P^s'>^>^) withp(.g)>0, "^geG, (135) 

is a sufficient condition for the existence of a mass gap. Indeed if we define p[G] = min p(g) and 

geG 

take the limit => oo, the r.h.s. of poip converges to 0(1) x c'^I'^l". thus p[G] gives a lower 
bound for the mass gap. It is also obvious that p[G] < follows \i H <Z G. {0{g)\V]) would 
be calculable by, for example, Monte Carlo Simulations and weak-coupling-, strong-coupling-, 
1/d- and 1/iV-expansions. 

• We proved the theorem on a square lattice, but it can be easily generalized to other lattices 
such as a triangular lattice. 

• Note that theorem [5] is derived with no knowledge of the interaction except for the reflection 
positivity. The strength of interaction can be made anisotropic, since it respects reflection 
positivity. And it is also correct in arbitrary dimensions. 



3.3 Examples 

There are many classes of lattice systems to which our theorem is applicable. One is the class of coset 
models^ where the field takes values in G / H with G an arbitrary Lie group and H its closed subgroup. 
We present three explicit examples. 

Example 1. The 0{N) Heisenberg model. 

The partition function is given by 

Z = f l[d$y exp(^(3Y.$--^-+i^)' $e S""-' =0{N)/0{N -I). (136) 

The most natural definition of a correlation function is 

r^(a;; n) = (j)^ ■ (j^x+nji- (137) 

In this case Z2 is obviously the smallest possible invariance group; of course, 0{N) itself can also be 
chosen. One can check that the requirements ([99|) . ()100p are met for both of them. If one is going 
to choose other arbitrary subgroup of 0{N), the correlation function should be redefined properly as 
explained below (|100p . Since the length of spin is normalized, the prefactor in poip becomes 1. 

Example 2. The CP^'^ model. 



The partition function is given b; 




^= / Il'^^i' ^^p(^IlTr[P,P,+^]), P eCP''-^ = U{N)l{U{l)xU{N -I)). (138) 

P is an N X N matrix and obeys P^ = Px, P^ = Px , Tr P^, = 1 . (An alternative way is to express 
Pij = z* Zj, where z € and = 1.) The action is invariant under the global transformation P^ 

IS See also ref.[34]. 
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fip^n^ with n e 5[/(iv)0 In two dimension, the model is beheved to possess a nonperturbatively 
generated mass gap for which, however, no rigorous result is available. Defining an appropriate 
correlation function needs some care; we define 



r^(x; n) = Tr ^ [P, - ^][P.+nf. - ^]| = Tr [P.P,+„,,] - j^, (139) 



where In denotes the unit matrix of size N x N. Then it is easy to confirm ()100|1 for G = SU{N) 



SU{N) 



N \^l JsU(N) 'J N ^ ^ 

^ MfeL_i!^^0. (141) 
N N ^ ' 



Example 3. G x G principal chiral model (PCM). 

Here G is an arbitrary compact group. The partition function is given by 

Z ^ j X{dUy c^vip ^°Tr [C/,(C/,+^)-i] , C/, e G. (142) 

yeA y x.fiCA J 

This model is invariant under a global G xG transformation U ghUg^^ . The corresponding twisted 
partition function is given, in agreement with the original definition p2p . by (Vtwall) 

Z{9) = I WdUy e^vip I^*^Tr [C/,(C/,+^)-i] + /? ^ RcTr \U,{gU,+ -^)-^] j , geG. 

veA \ x.p,cA\V x^iCV J 

(143) 

The correlation function is defined as 

r^{x;n)=XR{UxU-^„p, (144) 

where R denotes an irreducible unitary representation of G. If we choose G itself, then (jlOOp is satisfied 
if and only if i? is a nontrivial representation. This is a direct consequence of the Schur orthogonality 
relation in representation theory; see Corollary 4.10 in ref. [33]. If G = SU{N) and we choose Z^, 
then poop is satisfied iff R has a nonzero iV-ality. Anyway the prcfactor in (|10ip becomes 1. We 
emphasize that G need not have a nontrivial center. 

As an example other than cosct models, we only present 

Example 4. SU{Nf) x SU{Nf) linear sigma model. 

Usually this model is used to study the chiral phase transition [35]. On the lattice, the partition 
function is given by 

Z= I \{d% exp(/3i^ReTr(<i>,$t^^)-/32^Tr($t$,)_Ai^ [Tr(<I>t$,)]'-A2 5]Tr(<i>t$, 

ySA x,^ x^K xGA sGA 

(145) 

GM(iV/,C). (146) 
This model is invariant under ^ c*"gi$gjj with c*" e Ua{^) and gL,gR G SU{Nf). 



''The true symmetry is SU(N)/Zj^ , since SI € Zm does not change Px at all. 
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It might be the case that whether the upper bound il01\) gives an exponential decay or not depends 
on the choice of G, a point worth further study. 
* * * * 

Although we proved theorem [5] only in the case of the nearest-neighbor interaction, we can 
prove it even in the presence of non-nearest-neighbor and multi-site interactions if some 
appropriate conditions are satisfied. To make the argument concrete, let us consider the SU{N) x 
SU{N) PCM and suppose that Zjy C SU{N) was chosen as a symmetry group for theorem [S] There 
are two crucial conditions: one is site-reflection positivity and the other is that twists and their algebra 
( see 13. 1 p should remain well-defined. Here is a partial list of possible extensions (on a square lattice): 

1. The multi-site interaction term between four variables on the same plaquette, 

RexB.{U,ul^-U,^+f.+oUl^,), (147) 
can be added to the action without spoiling theorem [5] if the 7V-ality of i? is 0. 

2. The non-nearest-neighbor interaction term Re XR{UxUl_^2fj.) can be added to the action if the 
A^-ality of i? is and the coefficient in front of it is positive. (If the distance is larger than two 
lattice spacings, or if the coefficient is negative, then the site-reflection positivity becomes hard 
to prove.) 

3. The non-nearest-neighbor interaction term Kc xr{UxUI_^p^_^_0), fJ, ^ i', can be added to the action 
if the iV-ality of i? is 0. 



Note that our proved inequality may fail to be useful in phases other than the disordered phase, 
even though it is correct in any phases. Consider a spin system with a global symmetry group G in 
a c?-dimensional box of size Li x ■ ■ ■ x Ld whose boundary condition is twisted by g G G in the x^- 

d 

direction and otherwise periodic. Let denote the twisted partition function and set L± = JJ^ Lk- 

k=1 

Based on our experience in gauge theories, we generally expect following behaviors of jZ in various 
phases: 

Z^ jZ « cxp(— .tLj^ exp(— j/Li)) (Disordered phase) 

Zs/Z w exp(-zL_L) (Ordered phase) (148) 

Z^ jZ w exp(— wLj^/Li) (Massless phase) 

where x, y, z and w are functions of the coupling constants and the choice of g. Letting a denote the 
value of (1 — Z^ jZ^I^^ in the thermodynamic limit, we have a ~ c^^ in the disordered phase, while 
a = 1 in the other cases so that the r.h.s. of (jlOip tends to a constant independent of n. So the 
bottom line is that algebraic decay of correlation function cannot be inferred from the behavior of the 
r.h.s. in general. 

Finally we comment on the formal difference between the inequality derived by Kovacs and 
Tomboulis in refs.[l9l|20] and ours in the two-dimensional SU{2) x SU{2) PCM. Their result is 

/p / \v ^ Zn{-\-, + ) — Z„(— , — ) (-[AOl^ 

[L,[x; n)) ^^^^ _ ^^^^^ ^ ^^^^^ _^ ^ Z„(-, +) + Z„(-, -) ' ^'^^^ 

where Z„(Ti,r2) (ri.2 = ±1) is the partition function on the lattice of size nx n with a twist 
operated in the x^-direction. 
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On the other hand, our resuh f poip with G ~ Z2) gives 
(r^(x;n))<2{l-(^o[V])r/^" 

ZlS+.+)-ZlS+,-) 



(150) 
(151) 



Note that in the latter, both r.h.s. and Lh.s. are estimated on the lattice of size L^. In both formulas 
the correlation function is assumed to be in the fundamental representation. Although they look 
different, both relates the exponential suppression of the wall free energy (in the thermodynamic 
limit) to the mass gap, thus their physical contents are totally consistent. 

3.4 Demonstration in ID PCM and 2D square Ising model 

Let us verify the proved inequality (jlOip explicitly in the G x G PCM in one-dimension, with G an 
arbitrary compact group. The partition function of the model on a periodic chain of length L is given 

by 



» L L 

Za= Y[ dUk exp (/3 Tr (C/,C/^\)) , /? > 0, C/, G G, 

k=l i=l 

„ L L 

^ HdUk l[[Y,drFrXriU,U-^\)\, 

J I 1 A 1 ^ 



(152) 
(153) 



where runs over all irreducible unitary representations of G. Er = F- > follows from the 

r 

reflection positivity and reality of the action. Let us define = Fr/Fo for later convenience (0 
denotes the trivial representation). Straightforward calculation yields 



Z^=Y,dl{Fr)'^. 

r 

The twisted partition function is similarly given by 

. L L 
Zi= \{dUk [Y.d,.Fr^Xr'{9UlU^^)\\{[Y,drFrXr{U^U-^^) 

k=l r' 1=2 r 

= Y.'^rXr{g){Fr)'' , 



(154) 



9 G G', (155) 
(156) 



where G' C G is an arbitrary subgroup of G. (Z^ reduces to Za for g = 1, as it should be.) Hence 
we get 



hm {1- / dg{0{g))] 



n/L 



lim 

L — >OQ 



{Cr'Y 



Y^dliCrY 

r 



(157) 



Here is defined as a sum over all representations of G which 

r 

defined as the largest one among {c,- | r is nontrivial w.r.t. G'}. 



(158) 

are nontrivial w.r.t. G' and Cr' is 
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Next we define the correlation function as r(n) = Xr{UqU,^ ^). Then (|100p requires R to be 
nontrivial w.r.t. G. After straightforward calculation we get 

inn)) 1 . 



(r(o)> dn 



iXRiUoU-')) = (159) 



in the thermodynamic limit (L — > oo). Since c^j; < c^' is obvious from their definitions, we conclude 
from (|158|1 and (|159p that the inequality (jlOip certainly holds at least in the limit L — > oo. El 
* * * * 

As a next example let us take the two-dimensional Ising model on a square lattice. The partition 
function of the model is given by 

■^A = / duk exp ^ ^(00-^^,0-^+1,^ + 6o-^,.o-^.,^+i) , (160) 

where a^^y is the Ising spin located at the site (/x, v) and j da = — . Periodic boundary conditions 

are imposed so that ai^y = (tli+i.u, 0-^,1 ~ (^h,Lq,+i- We assume a > 0, 6 > 0. Let us focus on the 
high temperature (disorder) phase of the model. 

The exact asymptotic form of the two-point correlation function is known [36] and the mass gap 
(or inverse correlation length) M = 2(a — 6), where a is the dual temperature defined by 

sinh2a sinh2a = 1. (161) 

b is defined in the same way. 

To estimate the free energy of walls we need explicit formulae for twisted and untwisted partition 
functions. Here we use the expressions due to Kastening [37] , which in out notation read 

Zj^ = l[2sinh(2a)]^i^^/2x 

^1 M , TT To ■ 1 ^1 M , TT To 1/^1 M TT To • 1 ^1 



n [2cosh(i^72.-i)] +n [2Binh(^72.-i)] +n [2-Bh(^72.-2)] -n [2-nh(^^^^^ 

Kk=l k=l k=l k=l 

(162) 

Zl^^ = y[2sinh(2a)]^i^^/2x 
'ft [2cosh(^72.-i)] +ft [2sinh(^72.-i)] -ft [2cosh(^72.-.)] +ft [2sinh(:^ 



, fc=i fc=i fe=i fc=i 



72fc 

(163) 



•* is the twisted partition function; more precisely, it is a partition function on a lattice which 
is antiperiodic in x^-direction and periodic in a;^-direction. (Note that this boundary condition is 
equivalent to the existence of a wall wrapping around a periodic lattice in x^-direction.) 7^ > is 
defined by 

cosh 7fc = cosh 2a cosh 26 — cos— — sinh2a sinh26. (164) 

L2 



^® mOll l should hold for finite L too, but expressions of both sides of IHOH become highly eomplicatcd for finite L and 
verification seems to be hard. 
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The inequality to be checked, namely poip for the square Ising model, is given by 



1 / Z 



n/Li 



(aoa„)A<2|^ l^l-^jj . (165) 

(. . . )a denotes the expectation value measured on a finite lattice {— A). Similarly (. . . )oo denotes an 

expectation value in the thermodynamic limit. For simplicity we calculate not — — — but 

Za 



1 - J]^ tanh i — !-72A;-2 



7. y^-'^ ^2 



-72/c-l 



^'^^^^ cosh (-^72fc_i) L2 sinh(^ 

— f^n^^^ — 

fc-l cosh (^^72fc-2j cosh (^-^72fe-2 



(166) 



Considering thaim < 70 = 2(a — 6) is the smallest among {7*:}, we obtain, after some algebra, 

.l^^o.fl-¥) ' ^=exp{-(7o + f eVi)-V)}. (167) 



fe=0 



2L2-1 , 
Since E (-l)''+'7fc = o(— ) for L2 > 1, we get 



/c=0 



lim lim 

hi — >00 h\ — >OC> 



h.s. of (UnSl) = 2e"^°". (168) 



We compare this result with the asymptotic form of the exact two point function in the high temper- 
ature phase [55] : 



1 



1 



(to(t„)oo = /(a,6)^c-^«" X l + 0(— ) for n»l, (169) 



n 



where the factor /(a, h) is independent of n. 

(jl68|) and (|169p tells that the exponential decay rates of both sides coincide exactly for all values 
of (a, 5) when the system is in the disorder phase. This result suggests that our inequality might be 
a rather accurate one in general. 

3.5 Demonstration in 2D triangular Ising model 

Our next example is the two-dimensional Ising model on a triangular lattice. In the square Ising 
model, exact expression for mass gap was already known, thus the value of our theorem is obscured. 
However, as for the triangular Ising model, an exact expression for mass gap is not known except for 
special cases, so (unlike in the previous section) the results we give in this section are essentially new. 

We consider a lattice of size Li x L2 with periodic boundary conditions!^ See fig[TO]for an example 
with Li = 6 and L2 = 3; upper and lower edges painted in blue should be identified, and also right 
and left edges painted in red should be identified. (This lattice is the one that appeared in the seminal 
work of Houtappel |38| in which an analytic formula for the triangular Ising model was obtained for 



^^0 < 70 stems from the fact that now the system is in the high temperature (disorder) phase. 
^"Li denotes the number of triangles. It is not the actual length of the lattice. 



26 



the first time.) Another triangular lattice commonly used in the literature is depicted in fig llll where 
edges are again colored for the purpose of indicating the periodic structure of the lattice. It is easy to 
prove that these lattices are equivalent if and only if Li is a multiple of 2L2- See fig[T2]for illustration 
of this fact. Numbers arc written to guide the eye; edges assigned with the same number should be 
identified. Hereafter we will assume this condition, but this is only a technical assumption and not 
essential for we will be interested in the limit Li oo. The reason we did not start with the lattice 
in fig llll is because it does not allow for simple use of reflection positivity. 




Figure 10: A triangular lattice of size Li x L2 
with periodic boundary conditions; the red and 
the blue ends are identified, respectively. This 
lattice is symmetric w.r.t. each of vertical axes. 




L' 



Figure 11: A triangular lattice constructed 
from a square lattice by addition of diagonal 
edges. Its periodic structure is indicated by 
coloring as in figfTUl 




Figure 12: An illustration of the fact that those 
lattices given in figs [10] and [11] are equivalent 
iff Li is a multiple of 2X2- Here Li = 6 and 
i.2 = 3. 




Figure 13: A Z2-twist represented by a blue 
loop on the dual (hexagonal) lattice. On the 
original (triangular) lattice, it is represented 
by a stacked set of links (colored in red) with 
couplings of opposite sign. 



Let us define a twist on a planar triangular lattice. A twist is a closed loop on the dual lattice, and 
the dual of a triangular lattice is a honeycomb (or hexagonal) lattice as shown in figdJl Note that 
the blue line in fig[T3]is a closed loop owing to the periodic structure of the lattice. It becomes, on the 
original lattice, a stacked set of links with a coupling constant of opposite sign, which is depicted as 
a set of red links in fig ll3l Note that introducing a twist to a periodic lattice as in fig ll3l is equivalent 
to imposing an anti-periodic boundary condition in the horizontal direction. 
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Let us remember that it is not the number of the waUs but rather the number mod 2 of them that 
is physically relevant. This 'Z2 conservation' of walls is a direct consequence of tr^ = 1 in the present 
model, and we can show it explicitly by a sequence of changes of variables a —a. For instance, the 
partition function containing one wall and that containing three walls agree completely as illustrated 
in fig ll4l in which red segments represent twisted links. 




Figure 14: The partition function does not differ for any odd number of twists, owing to the Z2 
conservation of the twist. 

Exact expressions for partition functions of a planar triangular Ising model with various boundary 
conditions were derived by Wu and Hu via 'Grassmann path integral method' [35]. Let Za {Z^j^ ■*) 
denote the partition function with periodic boundary condition in both directions (with periodic in 
vertical and anti-periodic in horizontal direction), respectively. On a triangular lattice, three different 
couplings can be defined in each directions, so let us introduce Ji as the coupling constant on vertical 
bonds in figs llOllf] and J2,J3 the other two. Reflection positivity however requires J2 = J3 (= J). 
Introduce 

ii EE tanh(Ji/kBr), i = tanh(J/kBT). (170) 

Our convention is such that > corresponds to ferromagnetic coupling. On the (t, ti)-plane, there 
is a line which corresponds to T = Tc and we will call it the "critical line" in the following. Under 
the change of notation Li ^ N and L-z M, the result due to Hu and Wu for this case reads 



ZA = l[2cosh3(/3J)]''"' 
4-^ = y[2cosh3(;3J)]^^" 



Qi 1 

2 ' 2 



+ Qi Q + fig, 1 - sgn (T - Tc)QoM 



2 ' 2 



■Qi Q - % 1 + sgn (T - Tc)Qo,o 



(171) 
(172) 



where 



p=0 q=0 



2^(p + /i) 
M 



-Acos^"(^.+ ")-Acos^^"(^ + ^) 2.(, + .) 



Ao = {l + tHiY + (ti + ty + 2t\l + h) 



A 



N 



M 



N 



B 



2\2 



with Tc the phase-transition temperature. Using the formulae given above, we can show 
Lemma 6. In the disordered phase { T > Tc ) we have 

lim hm 1 - = e-P, 



1/2 



(173) 
(174) 



(175) 
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where 



p ^ cosh-i (^^) > 0, (176) 



The order of two limits in ^75^ must not be changed. 



--<.<1 



In the above, > 1 and \B\ < 1 arc implicitly assumed; these can be shown for every 



(t, ti) E (—1, 1)'^ by elementary methods. (Note that ^y-rr- = 1 defines the critical line.) The proof of 

theorem ini is elementary but technically cumbersome, which we relegate to the appendix. 

To gain an intuitive understanding of the above result, let us see fig llSl in which the projection of 
1 — onto the (ti,t)-plane is drawn. The black region corresponds to the (anti-)ferromagnetically 
ordered phase, while the colored region to the disordered phase. Brighter color represents larger p, 
hence larger mass gap. Fig [15] clearly shows a symmetry under t —t; this is a manifestation of the 




-1 -0.5 0.5 1 

Figure 15: The projection of 1 — c^p onto the (ti, i)-plane. The black region corresponds to the (anti- 
)ferromagnetically ordered phase and the colored region to the disordered phase. Brighter (darker) 
color represents larger (smaller) p and especially p diverges on the t = line. The boundary of the 
colored region signifies the critical line. 

well-known fact that the triangular Ising model is invariant under simultaneous sign reversal of any 
two of Ji , J2 , J3 . 
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The inequality of our primary interest, namely (jlOip for the triangular Ising model, read; 



21 



(ao.„U<2|l(l-^)|' . (178) 

Letting il/ — » oo after — > oo, we obtain 
Theorem 6. 

(aoa„)oo < 26-"". (179) 
The above is the main result in this subsection; p is a rigorous lower bound of the true mass gap. 
We should keep in mind that the l.h.s. of p78p is a pair correlation between two spins on the same 
horizontal level as depicted in fig II 01 the two spins are not on the same lattice axis. 

Since the exponential decay rates are equal for both sides of the inequality in the the square Ising 
model, it is natural to expect so in the triangular Ising model too. The asymptotic correlation between 
two spins on the same lattice axis in the triangular Ising model was derived by Stephenson for both 
ferromagnetic and antifcrromagnetic couplings [40j . However, the asymptotic correlation between two 
spins off the axis is not found in the literature. We conjecture as follows: 



Conjecture. 

p is equal to the true off-axis mass gap for every (i,ti) in the disordered phase. 



The most straightforward way to test the conjecture would be to measure the mass gap directly 
via Monte Carlo simulation. However, as already seen from (|169p the exponential falloff generically 
receives power law corrections (the so-called 'Ornstcin-Zcrnike' decay [41j ) which makes a reliable 
fitting difficult. To evade this hamper would call for sophisticated methods such as the Monte Carlo 
Transfer Matrix Method |42]. A numerical check of the conjecture therefore seems to be a highly 
nontrivial task, and we defer it to future work. 

In a special case, analytical test is possible: when ti = the model reduces to the isotropic square 
Ising model and the off-axis correlation function reduces to the diagonal correlation function. From 
(|T751) it follows that 

Pk=o = cosh-' Am-t') \t\<tc^V^- 1, (180) 

which completely agrees with the exact diagonal mass gap obtained by Cheng and Wu in 1967 [36j . 

Further insight is gained by considering the isotropic case t = ti. Since the two-point correlation 
function in this case is expected to be approximately isotropic (except for sign in antifcrromagnetic 
case), it seems reasonable to compare e"'' with e"*^^/^^™ where m is the exact on-axis mass gap |1D]. 
(■\/3/2 is a geometric correction factor.) 

Fig [TBI depicts the graphs of e"'' and Q-i^f^)™ against t e [—1, 2 - -^3]. For t > they agree 
quite well; their nonzero difference is hardly discernible to the eye. For t < agreement is still not 
bad. To say the least, the comparison suggests that p be fairly close to the true off-axis mass gap and 
supports, rather than defies, the conjecture. 

It is readily seen from (|177p . (|176p that c~p is a nonanalytic function in the region of negative 
ti, which, assuming the validity of our conjecture, implies non-analyticity of the mass gap. Such 
an exotic possibility definitely deserves further study. Strictly speaking, however, there are different 
possibilities that cannot be denied here: for example it could be the case that p equals the true mass 

Actually, we originally proved 11011 1 on a square lattice, but the whole procedure of the proof goes over to the case 
of a triangular lattice almost unchanged. 
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t 

Figure 16: c~'' and c~'^/^)™ arc plotted against t in the case of isotropic case {t = ti). Their 
agreement is remarkable especially at t positive. 

gap only when ti > 0. In the latter case, non-analyticity of p does not signify that of the true mass 
gap. 

Let us end this subsection by invoking the effectiveness of our approach. Although the circumstance 
concerning our conjecture is rather moot, it can be safely said that our result in this subsection is 
essentially new to the extent that it rigorously gives a lower bound for the still-unknown off-axis mass 
gap of the triangular (both ferromagnetic and antiferromagnetic, both isotropic and anisotropic) Ising 
model in the disordered phase. 

3.6 Strong coupling analysis 

In this section we show, using the convergent strong-coupling (taken as synonymous with high temper- 
ature) expansion, that both sides of the proved inequality poip have an identical exponential decay 
rate at long distance as long as the on-axis correlation function is considered!^ The proof is valid in 
any dimension and makes no use of reflection positivity. Since the corresponding result in LGT has 
already been derived by Miinster [2S] (as mentioned in section and since no essential difficulty 
arises in extending his proof to the case of spin models, we shall be brief here and only try to sketch 
the main idea behind the approach. Implications of this result to our conjecture will be discussed 
later. 

For simplicity of exposition let us consider the isotropic square Ising model and its on-axis corre- 
lation function (though our argument is readily extendable to more general models such as PCM). A 
precise statement of the claim goes follows: as long as the size of the lattice is larger enough than n, 
the strong coupling expansion (SCE) of p and m are identical at least up to order n. (Our notation is 
such that the definition of p is in (|175p . m is the mass gap, are the same as in section 

[SHand t = tanh(J/kBT).) 

High-temperature behavior of correlation functions in Ising-Uke models have been studied by many authors in a 
variety of methods; see ref. 1431 . for example. It is worthwhile to note that a majority of existing studies deal with neither 
the off-axis correlation function nor the case of an antiferromagnetic coupling. Hopefully a partial understanding of 
this fact will be gained through the discusions in this subsection. 
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Let us begin with the expression Z = + t<^i<7j)- Expanding Z into sums of diseonnected 

{t} i^] 

loops and then taking the logarithm, we have logZ = ^^i'''' with 7 any connected loop and I7I the 

7 

perimeter of 7. Using similar expression for Z^~^ we obtain log — ^— = —2 t''''' where 5 is the set 

of loops wrapping around the lattice in x^-direction for odd number of times. Since we are interested 
in the limit \t\ <^ 1. it is sufficient to consider only such loops that wind around the lattice in x^- 
direction only once. Factorizing the degeneracy factor due to translational symmetry in x^-direction, 

we have log ^ = —2^2 ' definition of S' should be obvious. 

76S' 

It is clear that the leading contribution, of order 0{t^^), comes from a straight line extending 
in .T^ -direction while the subleading contributions come from loops which are formed via addition of 
some 'decorations' to the leading line. Dividing by the leading contribution and taking the logarithm 
will single out contributions of connected decorations, which is proportional to Li owing to the 
translational invariance of the straight line. Thus we find exactly the behavior (|148p in section [J 



log 



ocLi. (181) 



On the other hand, the on-axis two-point correlation function {a'xO'x+r)oo can be written as a sum over 
contributions of lines connecting tri to (T2 , whose leading term comes from a straight line extending 
between ai and (72 and subleading terms from zig-zag lines that descend from the leading one through 
addition of decorations. In this way we see that the SCE of 

hm log[f-Llog^) A^iI/Li {=-p-\ogt) (182) 



L\,L2 — ^00 L \ L2 Z 

is identical , term by term, to the SCE of lim log \ {<Tx<Jx+r) 00 /t''] / {— — logt). Hence m ~ p. 

r — >oo 

The argument above is valid for various other models as long as on-axis correlation functions are 
concerned. Then it is natural to ask about off-axis correlation functions. (This is the case relevant 
for the conjecture.) From fig[TO]it is easily understood that the leading contribution to the SCE of 

— log does not come from a single straight line: instead it comes from ( ^ \ different loops, 

L2 Z \Li/2J 

all of the same length Li. So we now have a number of different ways to see a given higher-order 

loop as a sum of any one of the leading- order loops and a decoration added to it! This implies that 

1 z^-^ 

the counting of diagrams appearing in SCE of — log — — — (and of off-axis correlation function, too) 

L2 Z 

is immensely complicated. We even face another problem: since most of the leading-order loops have 
no translational symmetry, it becomes a nontrivial task to show (|18ip . For these reasons we cannot 
give a mathematically rigorous proof of the conjecture even at sufficiently high temperature. 
Some caveats are in order. 

• First, remember that we did confirm m = p for the diagonal correlation function in the square 
Ising model ( (|180p and the accompanying discussion). This fact implies that our inability to 
prove (the very existence of m, p and) the equality m = p for off-axis correlation function in 
SCE approach does not itself constitute a disproof. 

• One may be tempted to argue that, since an exact one-to-one correspondence between the 

1 z^-^ 

diagrams for SCE of (cricr2)oo and those for SCE of — log — -— exists, m = p would readily 

L2 Z 
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follow even in the ofF-axis case if we presume the existences of both m and p. This reasoning is 
however incorrect, because the exact correspondence is present only in the on-axis case. (This 

.1 z^-^ 

point is quite nontrivial.) An example of a loop that appears in SCE of — log — ^i— as one 

L2 Z 

of the leading contributions but has no counterpart among the diagrams in SCE of (cri(T2}oo is 
shown in fig |17l (Although somewhat counterintuitive, this loop has a minimum perimeter to 
wind around the lattice in horizontal direction.) 

Since we have to send Li — > 00 before L2 ^ 00 in estimating p, those loops can never be 
neglected. From this point of view, it is rather surprising that m = p holds in the ofF-axis case 
of the square Ising model. 

• It seems worthwhile to note a qualitative difference between ferromagnetic and antiferromagnetic 

cases. The so-called random path representation ((Ti(T2) = i'''' , where the sum runs over every 

7 

path connecting ai to a2, loses its probabilistic interpretation when t^'*^ < 0. This actually 
happens in a triangular Ising model with antiferromagnetic couplings. In such a case it is 
impossible to apply fertile probability-theoretical methods, so that a particular difficulty is 
envisaged in settling the conjecture when the coupling is antiferromagnetic. 

Summarizing above, we do not have a definitive answer as to the validity of the conjecture even at 
sufficiently high temperature, and further analysis on this topic will be hopefully reported elsewhere. 



Figure 17: Drawn in red is an example of a loop wrapped around the lattice in vertical as well as in 
horizontal direction. Such loops appear numerously if the horizontal size of the lattice is much larger 
than its vertical size. 



4 Conclusion 

In this paper we generalized the inequality of Tomboulis-Yaffe in SU{2) LGT to SU (N) LGT and also 
to general classical spin systems, together with a detailed analysis of basic properties of non-Abclian 
twists. Our result is obtained essentially on a finite lattice and gives a rigorous upper bound of a 
Wilson loop and a two-point correlation function. An intriguing point is that the inequality obtained 
for spin models does not require the center of the symmetry group, so they can be applied e.g. to 
G X G PCM with centerless G. This point seems to be a progress compared with preceding studies 
in which the center was perceived as special without physically convincing motivatiorj^. 

^^Some arguments that allege the speciality of the center do exist in the literature, but they seem to be rather subtle, 
as discussed attentively in ref. |21|. 
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Furthermore we demonstrated our result explicitly in some solvable models and found in the 
square Ising model that the obtained lower bound of the mass gap is equal to the exact one. We also 
calculated the off-axis mass gap in the triangular Ising model for various couplings, but this time the 
exact mass gap is not known and direct comparison is impossible. Wc conjectured that the bound is 
indeed saturated and pointed out that the conjecture implies the non-analyticity of the mass gap. We 
have tested its validity in several ways, including strong coupling analysis, but a definitive conclusion 
is still lacking and is left for future work. 

At present the mechanisms of the quark confinement in non-Abelian gauge theories and the mass 
gap generation in non-Abelian spin models still remain elusive, and we hope that our result will be 
useful for further clarification of the issue. 

Note added 

After this work was completed, we learned that C. Borgs and E. Seller had already obtained a 
result very similar to theorem [2] of this paper; see Lemma II. 8 and the accompanying discussion in 
ref. [440. But since it links the Polyakov loop correlator and not the Wilson loop with the electric 
flux free energy, it is not quite the same as ours. It does, however, already imply the 't Hooft's string 
tension is less than or equal to Wilson's (sec (11.48) and (II. 50) of ref.jH]). Finally wc note that their 
results hardly overlap with ours in section [31 
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Appendix A Proof of Lemma [6] 



For T > Tc, ([T7T|) and yield 



f2i 1 -\- r2i Q 



Let us define 0^{^l,p, M) > by 

cosh6'f (^,p,M) = 



1 — B cos 



M 



2 A cos 



tt{p + fi) 



M 



(183) 



(184) 



It is tedious but straightforward to show that the minimum of coshfl^ as a function of —1 < 
cos ^^^^^^ < 1 is given by (see ([rfT]) ). and that > 1 for every (t,ti) e (0, l)^ with 

equality on the critical line. After elementary calculations, we find 



JV 



p=0 



M 



(185) 



^'^ I thank E. Seller for kindly pointing out this fact to me. 
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for 1. Since (|185p has no dependence on i/, 

n 



lim ..... 

N^oo \li 1 N—>oo iZn i 

2 ' 2 



(186) 



Next, using pSSp we get 

i±.o.p{iVi:,-l)««/,= (^)} t„r«»l, 



r^i 1 

2 



2A/-1 



A;=0 



with 



/f (x) = y fog|l -Bcos(27rx) + (l - Bcos(27ra-))^ - (2Acos(7rx))^ |. 



(187) 



(188) 



^0 

Let us investigate how fast — — — converges to 1. Using (|185p we can show 



n {l-4exp(-iV0f)} for TV » 1. 



(189) 



p=0 



Define 6* as the smallest of {6*^(0, p, Af)}p. Then (fT^ simplifies to 



'0,0 



^0 i 



1-Ke 



-NS 



for TV > 1. 



(190) 



ii' is an integer G {4,8,12,16}, dependent on A, B and M. Substitution of (fTSS)) . (fW|) and ([TOO)) 
into (fTB5)) yields 



1 



(-) _ 2M-1 

2 exp I - iV 



Since 9 = cosh 



m 

\A\ 



O 



0+ E 

A:=0 
2A/-1 

L) and Y: (-D^+Vi 



fc=0 



2A/ 

k 
2AF 



)]} 



for TV > 1. 



(191) 



= O 



1 



for M > 1, wc find 



lim lim 1 ^ 



1/N 



p = cosh 



1 (9{B)_ 
\A\ 



(192) 
(193) 



which is the desired result. 



□ 
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